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We present a quantum field theoretical approach based on the Lehmann-Symanzik-Zimmermann 
reduction for the multi-photon scattering process in a nano-architecture consisting of the coupled- 
resonator arrays (CRA), which are also coupled to some artificial atoms as a controlling quantum 
node. By making use of this approach, we find the bound states of single photon for an elemen- 
tary unit, the T-type CRA, and explicitly obtain its multi-photon scattering S-matrix in various 
situations. We also use this method to calculate the multi-photon S-matrices for the more complex 
quantum network constructed with main T-type CRA's, such as a H-type CRA waveguide. 

PACS numbers: 03.65.Nk, 42.50.-p, 11.55.-m, 72.10.Fk 
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I. INTRODUCTION 

In order to realize all-optical quantum device&ii^, it is 
very crucial to explore the physical mechanism for the 
single photon generation, transport and one shot detec- 
tion etal. . Thus, we need a comprehensive understand- 
ing of the fundamental processes of coherent photonic 
scattering in the solid state based confined systems, such 
as the photonic crystal with artificial band gaps. Es- 
sentially, by coupling the system with an extra two-level 
system (TLS) to form a hybrid system, for controllably 
transport of photons the basic element is a quantum node 
or quantum switch'^. It was abstracted as the so-called 
photon transistor most recentlyi*^. For quantum infor- 
mation processing, such quantum node can coherently 
control the quantum state transfer in some quantum nct- 
work^^^. Actually, to manipulate the coherent transport 
of photons, the quantum node is tunable so that it can 
behave either as a perfect mirror totally refiecting pho- 
tons, or as an ideal transparent medium allowing photons 
to pass throughly. Theoretically, the quantum node for 
single-photon in all-optical architectures was extensively 
studied in one dimensional waveguide by making use of 
the standard scattering approaches^iSrJ^. 

The photonic quantum node is usually modelled as a 
localized TLS^*, which can be implemented as an arti- 
ficial atom, coupled to the photons transported in the 
coupled-resonator arrays (CRA)i^. The atomic parame- 
ters, e.g., the energy level spacing, are tunable to control 
the propagation of photons. Recently, based on this the- 
oretical model and its generalizations, the photon trans- 
port in the CRA systems has been studied for different 
purposes^i^ii^. Here, we would like to point out two re- 
markable issues: (a) if only one photon allowed to trans- 
port in the CRA with the artificial atom prepared at its 
ground state, the hybrid system can be described by a 
simple model with a single state coupled to a continuum, 
which is referred to the Anderson-Fano-Lee modeli^""— 
with single excitation; (b) the CRA can be regarded as 



the waveguide with nonlinear dispersion relation. It can 
be linearized in the high frequency regime in which pho- 
ton momentum k ^ ±T:/2a with lattice spacing a, to 
approximate the linear dispersion relation for the con- 
ventional waveguide. 

In the recent study^, the single photon transmission 
and reflection coefficients were calculated for the incident 
photon with any energy in the T-type structure (a quan- 
tum node coupled to a CRA, see Fig. 1), which demon- 
strate the novel lineshapes beyond the Breit-Wigner— 
and the Fano lineshapes^^ . This kind of investigation was 
carried out only for the case with single photon. Also, 
there are only quite a few researches on the two photon 
case, for which we mention an elegant theoretical ap- 
proach for two photon scattering^i^ based on the Bethe- 
ansata^"— . Moreover, we have to say that it sounds 
very difficult to prepare the system only with one pho- 
ton or two photons exactly, thus these previous studies 
need to be improved for multi-photon processes oriented 
by practical application. 

Actually, the study for multi-photon transport is very 
important to realize a practical all-optical devices. How- 
ever, these subtle approaches for single and two photon 
cases mentioned above^'^'i^, such as the discrete coor- 
dinate scattering theory and the Bethe-ansatz technique 
with fixed scatterers, are not feasibly generalized for the 
realistic multi-photon processes, even for two-photon or 
three-photon processes. Therefore, our present systemat- 
ical approach based on quantum field theory is significant 
since it is obviously feasible and intrinsically natural for 
the generalization to multi-photon scattering processes. 

In this paper, we utilize the Lehmann-Symanzik- 
Zimmermann (LSZ) reduction^^ in quantum field theory 
to investigate the multi-photon transport in the complex 
CRA with some two- level scatterers. This method has 
been used to study the single electron inelastic scattering 
in Anderson model and Kondo modeP^'"^°. Here, we deal 
with the multi-photon scattering problem by studying 
the out-state of the scattered photons for an arbitrary 
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state of incident photons. In a middle stages, we cal- 
culate the multi-photon scattering matrix (S'-matrix) in 
details. With the diagrammatic analysis, we find that 
the basic element of the S'-matrix is a connected trans- 
fer matrix (T-matrix), which can be obtained from the 
well-known LSZ reduction formula about the photonic 
Green's function. From the explicitly achieved expres- 
sions of photonic out-states, we analysis in details quan- 
tum statistical characters of photon transmission in the 
situation with many photons. We find that, in the tight 
binding CRA of T-type, there exist the single photon 
bound states. As a test, two photon transport in the T- 
type waveguide is re-considered, and our obtained results 
accord with the recent woks^ using the Bethe-ansatz, 
which verifies the results based on the LSZ reduction ap- 
proach are valid. As the development, the three photon 
scattering is studied, and the outgoing states of the three 
photon are given by this approach. Our present investi- 
gation mainly based on these results, can be regarded as 
a substantial development for its particular emphasis on 
the multi-photon scattering. 

In practice, the T-type photonic element we mentioned 
above is the basic block to constitute a complex quantum 
network coherently transferring photons in a controllable 
fashion. A slightly complicated illustration of such archi- 
tecture is the CRA waveguide of H-type. In this paper, 
we also study two photon scattering processes in the EE- 
type in details. 

The paper is organized as follows: in Sec. II and III, we 
model our hybrid system for multi-photon transport and 
present the scattering matrix based on the LSZ reduction 
approach; in Sec. IV, we show that there exist the single 
photon bound states in the tight binding T-type CRA; 
in Sec. V, we study the multi-photon transport in the 
T-type waveguide; in Sec. VI, we discuss the two pho- 
ton transport in the H-type waveguide; in Sec. VII, the 
results are summarized with some remarks. 



II. SCATTERING MODEL FOR THE HYBRID 
SYSTEM 

A. Model setup 

In this subsection, we model the T-type CRA by the 
two-level atom coupled to photons inside CRA illustrated 
in Fig. [TJ The model Hamiltonian reads 

Ht = n |e)(e| -I- Yy^oa\a, - J{a\ai+i + H.c.)] 

i 

+vY,5,o{a\a- +h.c.), (1) 

i 

where the operator (t~ = \g) (e| denotes the flip from 
the atomic ground state \g) to the excited state |e) 
with the energy level spacing Q,. Here, J is the 
hopping constant characterizing the inter-cavity cou- 
pling in the tight-binding approximation; (a|) is the 




FIG. 1: (Color online) The schematic for the complex CRA: 
The T-type structure with CRA coupled to impurity (maybe a 
two-level scatterer)is the basic element to construct the more 
complicated architecture of quantum network. The red circles 
denote the two level impurity, while the blue dots denote the 
photonic coupled resonators. 

annihilation (creation) operator for the photonic single 
mode with eigen-frequency wq in the i-th cavity; V is the 
hybridization constant of the localized atom-photon in 
the 0-th site of CRA. 

In the momentum space (fc-space), the Hamiltonian ([1]) 
is re-expressed as 

Ht = ^ |e)(e| + y^^eka\ak 

k 

+ ;^E(4'^"+H.c.), (2) 

with the photonic dispersion relation 

Efe = Wo — 2Jcos fc , (3) 

in CRA of length L. Here, we choose the cavity spac- 
ing a = I. In the high energy limits k ~ ±7r/2 and 
= ttJ, the above dispersion relation is linearized as 
£k ~ which is the same as that in the conventional 
waveguide. Thus, we can use CRA to simulate the con- 
ventional waveguide in the high frequency limits. As the 
basic element of the complex quantum network as illus- 
trated in Fig. [TJ the two- level atom plays the role of the 
photon transistor to control the photon transmission. 

It is obvious that, when confined within the single ex- 
citation subspace, our model is the same as the models 
by Anderson, Fano and Lee. This observation motivates 
us to consider how to use various approaches developed 
previously for the models by Anderson, Fano and Lee, 
to deal with the coherent processes of our system, espe- 
cially with multi-photons. To this end, we will compare 
our model ^ with the Anderson, Fano and Lee model 
as follows. 
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B. Relation to Anderson model 

By neglecting the Coulomb interaction, the Anderson 
model is the same as the Fano model, so we discuss the 
similarities and differences between the model ([2]) and 
Anderson model. The corresponding Hamiltonian reads 

Ha^Y. ^kcl^^ka + Sdflja +Hv+Hu, (4) 
k,(T 

where Cka- (c^g.) is the annihilation (creation) operator of 
the conductive electron with dispersion relation Ek and 
spin a. fa- (/J) is the annihilation (creation) operator of 
the impurity /-electron with energy Sd and spin a. The 
Coulomb repulsive interaction of the impurity /-electrons 
is described by the Hubbard term Hu — Uf^fiflfi- The 
hybridization of a conductive electron with the localized 
/-electron is depicted by the mixing term 

Hv^vY^iclja + H-c). (5) 

Anderson even used this model ^ to investigate the 
localized magnetic state in the metal. By comparing the 
Hamiltonian (|3|) with the Hamiltonian we find that 
the atom and the photon in the model ^ can be con- 
sidered as the impurity /-electron and the conductive 
electron in the Anderson model, respectively. But there 
still exist some differences between the two models. First, 
the Hilbert space of the atom in the model ([2]) is spanned 
by two basis states, |e) and \g), but the Hilbert space of 
the impurity /-electron is spanned by four basis states, 
|0), /| |0), fl |0) and /j/| |0). Secondly, in the Anderson 
model the conductive electrons are fermions, while in the 
model photons are bosons. Thirdly, there is no Hub- 
bard term Hu in the model while in the Anderson 
model if there is only one /-electron the Coulomb interac- 
tion vanishes, i.e., the Hubbard interaction Hjj does not 
play any role. In conclusion, the two models are equiva- 
lent only in the single excitation case, i.e., one conductive 
electron with no /-electron or one impurity /-electron 
with no conductive electron in the Anderson model, and 
one photon with the atom prepared at the ground state 
or no photon with the atom prepared at the excited state 
in the model ([2]). Otherwise, in the multi-excitation case, 
the two models and ^ are obviously not equivalent. 

C. Relation to Lee model 

The Lee model^* with the Hamiltonian 

k 

describes a reaction process 

V^N + Ak, (7) 



with two heavy fermions (F, A'') of masses my and 
•niN- Here, the relativistic boson Ak with momentum 
k and rest mass /i possesses the dispersion relation ujk — 
yfc^+T?, and g is the three body coupling constant for 
the scattering of the three kinds of particles, V , N , and 
A. T. D. Lee used this exact solvable model to study 
the necessity of renormalization in quantum field theory 
even without the perturbation expansion. 

By comparing Lee model with the model ([2]), we find 
that the heavy fermions V and N can be regarded as the 
excited state and ground state of the atom, respectively. 
And the relativistic boson Ak can be considered as the 
photon in the model There also exist some differ- 
ences between the Lee model and the model ^ . Firstly, 
in Lee model, the renormalizations of the T^-fermion mass 
mv and coupling constant g should be taken into ac- 
count, and the Hamiltonian ([5]) is non-Hermitian. Thus 
we should investigate the unitarity of the S'-matrix care- 
fuUy^i^^. Secondly, the Hilbert space of heavy fermion is 
spanned by four basis states, i.e., |0), 4'v 1^)' 1^) ^^"^ 
V'yV'jv 1*^)- We conclude that the two models ([2|) and (|6]) 
are equivalent only in the single excitation case, i.e., one 
relativistic boson Ak with one heavy fermion N or one 
heavy fermion V with no relativistic boson Ak in Lee 
model, and one photon with the atom prepared at the 
ground state or no photon with the atom prepared at 
the excited state in the model 

In fact, it is turn out that the multi-particle scattering 
problems (such as N-69 scattering) in the Anderson and 
Lee models were both successfully studied by the LSZ 
reduction approach— i^i^^iM. In the following, we start 
with the model ^ to investigate the 5-matrix of scatter- 
ing photon in the hybrid CRA systems. For the further 
discussion, we first give a general formalism to study the 
multi-photon S'-matrix in the complex CRA by the LSZ 
reduction formalism. 



III. LEHMANN-SYMANZIK-ZIMMERMANN 
REDUCTION FOR PHOTON SCATTERING IN 
GENERAL 

In this section, we first briefly summarize the main 
results of the LSZ reduction approach in quantum field 
theory, and then we use them to give a general formula 
for multi-photon S-matrix in the T-type CRA. For in- 
vestigating the multi-photon scattering, in the first sub- 
section, we first define the n-photon S-matrix and 2n- 
point photonic Green's function, which are explicitly de- 
scribed by the Feynman diagrams (Fig. [5]), and utilized 
the functional integral to obtain the 2n-point photonic 
Green's function. In the second subsection, we use the 
LSZ reduction approach to obtain a general form of the 
n-photon S-matrix element by reducing the external legs 
(red curves in Fig. [5]) of the Green's function. 



A. Main results of the LSZ reduction approach 

Because the LSZ reduction approach relates to the 
photonic Green's function and ^-matrix, we need to give 
their definitions for discussing the LSZ reduction exphc- 
itly. The 2n-point photonic Green's function 



G 



= (e|TapXi'i)-«P„(4)«L(^i)-4„(^n)|0) : (8) 

is defined by the time ordering product T of the photon 
creation (annihilation) operator aj, (a^,) in the Heisenberg 
picture. Here, 6 denotes the ground state of the Hamil- 
tonian H of the system. 
The S'-matrix element^i^^ 



(/I 



{f\S\ 



(9) 



of photons is defined through the overlap of incoming 
state |z) and outgoing scattering state which are 

asymptotically free, and may contain multi-photon exci- 
tation, i.e., they are multi-photon states. In our case, the 
incoming state is 



|0) . 



In the present case, the S'-matrix is given by 



S = Texp[- 



+ 00 



(10) 



(11) 



with the atom-photon hybridization Hamiltonian 

H,nt = Viala- + H.c.) , (12) 

in the interaction picture. 

Using the diagrammatic analysis, we find that the 
Feynman diagram of S'-matrix element is constructed by 
summing up the contributions from all kinds of discon- 
nected diagrams as shown in Fig. Here, each of these 
disconnected diagram is made up of some connected di- 
agrams which describe the T-matrix elements. As the 
special cases, we consider the diagrammatic construction 
for the single photon and the two photon S'-matrix ele- 
ments by the connected T-matrix elements below. For 
the single photon case, the S-matrix element 



5p;fc 



5kp + iTj 



p;k ■ 



(13) 



is defined by the T-niatrix element for single photon, 
where k and p are the momenta of the incoming and 
outgoing photons. In this case, there exist two kinds of 
disconnected diagrams (a) and (b) as shown in Fig. [H 
For the two photon case, the S-matrix element 

^PlP2\kik2 ~ Sj92fc2 ~^ ^P2klSpik2 ~l" *-^PlP2 ;fcl ,^2 1 (1^) 

is reduced by the T-matrix element of two photons, where 
kr and p^. (r = 1, 2) are the momenta of the incoming and 
outgoing photons. In this case, there exist three kinds 




n : 





-|- other disconnected diagrams 



FIG. 2: (Color online) The Feynman diagrams for LSZ reduc- 
tion: The red line denotes the free photon propagator. The 
gray shade circles denote the S'-matrix elements and the blue 
shade circles denote the T-matrix elements. The Feynman 

cted 



(a) 



(b) 



FIG. 3: (Color online) The diagrammatic constructions of 
the single photon S-matrix element: There exist two kinds 
of disconnected diagrams (a) and (b). The red line denotes 
the free photon propagator. The gray shade circles denote 
the S-matrix elements and the blue shade circles denote the 
single photon T-matrix elements. 



of disconnected diagrams (a), (b) and (c) as shown in 
Fig. m Obviously, the multi-photon S-matrix is totally 
determined by the connected T-matrix, so we only need 
to find the photonic T-matrices. 
Fortunately, the intrinsic relation 



iT2n = G2nl[[2nG^\kr)G^\pr)] 

r=l 

between ?i-photon T-matrix element 

^271 ^pi^...^pn;kl,...,kn 1 

and photonic Green's function 

G2n = Gpj^^,,,^p^-ki,...,k„ I 



(15) 

(16) 
(17) 
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FIG. 4: (Color online) The diagrammatic constructions of 
the two photon S-matrix element: There exist three kinds of 
disconnected diagrams (a), (b) and (c). The red line denotes 
the free photon propagator. The gray shade circles denote 
the ^-matrix elements and the blue shade circles denote the 
two photon T-matrix elements. 



is given by the LSZ reduction formula, where 



Go(fcr) 



^0^ 



(18) 



is the Green's function of free photon and 
Gp-^....^p^-ki,...,kn is the Fourier transformation of 
Eq. ([5]). Here, the subscript os denotes the on shell limit 
w — » fffe. Finally, the multi-photon S'-matrix elements are 
determined by the photonic Green's functions entirely. 
In this paper, we use an elegant method, i.e., functional 
integrals, to establish the desired relation ([15]). This 
method has been used to solve many quantum impurity 
proble m^^i^" in the condensed matter physics. Here, we 
show how it works for the multi-photon transmission in 
this hybrid system. 



B. General formula for S'-matrix in the T-type 
CRA 



We utilize the generating functional Z to represent the 
full time ordering Green's function as 



pi,...,p„;ki 



5%^-^'n*pJvkt-Mk„ 



(19) 



where the generating functional 

exp{i[5 + fdt Y^ivlak + r/fc4)]} ,(20) 

k 

is defined by the action S = J dt L and the Lagrangian 
L = ifldtfe + iPgdtfg + I J2 aidtak - Ht • (21) 

k 

In Eq. (I^Ul) . we represent the two- level atom by the 
fermions fa- as 



|e) (e| - fife , 
\e){9\ = flfg, 



with the constraint 



(22) 



(23) 



This constraint arises from the fact that the physical 
space of the atom spanned by |e) and \g) is two dimen- 
sions, while the physical space of the fermions spanned 
by |0), ft |0), 4 |0) and |0) is four dimensions. By 
integrating the photon field and the fermionic fields f^ 
in Eq. (PO)) . the generating functional is obtained as 



dk 



-I I dujdk 
where the field variable is 

and the matrix 



oj - Sk + jO+ 
%-(w) 



27r cj - efe + iO+ 



(24) 



(25) 



is defined by the self-energy 



E(c^) = i?eE(c^)+i-p(^), 



(26) 



(27) 



of the atom and T — V^. Here, the real part of the 
self-energy is determined by the principal-value integral 
as 



„ „ / N I dk ^ r 
i?eE(c.) = / —P- 



2n £/£ — 



(28) 



and the imaginary part is proportional to the density of 
state (DOS) 



(29) 



6 



dk£k\ 



and Zi is the real root of the 



where Di{uj) 

equation e^- — uj. For the waveguide the self-energy 
is a constant. For the CRA the self-energy depends on 
the frequency uj and we have used Markov approximation 
to obtain Eq. (p4)) . 

With the helps of Eqs. ([T5|), dH]) and we can 

achieve the multi-photon ^-matrix by considering the 
Green's function G2„. In the following, we use the LSZ 
approach to study the multi-photon transport in the com- 
plex CRA. 



IV. BOUND STATES OF THE SINGLE 
PHOTON TRANSPORT IN THE T-TYPE CRA 



S- matrix element 
of the single photon 




The exact Green function 
of the two level atom 

(a) 



C0„ 



2J 

1 



b 



c(d) 
b 



(b) 



In this section, we consider the single photon trans- 
port in the T-type CRA. This problem has been consid- 
ered in Ref.^, but the photon bound states have not been 
taken into account. There are two interesting physical 
phenomena about the photon bound states, (a) Firstly, 
in the hybrid system the single photon bound states de- 
pict the states that the single photon is almost local- 
ized in the cavity containing the two-level atom. In this 
sense, the photon bound states can be used as the quan- 
tum information storage of single photon, (b) Secondly, 
as the basic element of single photon transistor, the T- 
type CRA is usually coupled with other CRAs through 
the two-level atom. These CRAs can be regarded as the 
scattering channels. It is proved that^^ the incident pho- 
ton in another CRA, whose energy is resonance on the 
bound state energy in the T-type CRA, will be totally 
reflected. This phenomenon can be considered as the 
one-dimension photonic Feshbach resonance, which can 
be used to simulate the Feshbach resonance in the atomic 
and condensed matter physics. The above reasons moti- 
vate us to study the photon bound states in the T-type 
CRA by using the LSZ approach. 

The T-type CRA is described by the Hamiltonian ©. 
Then Eqs. ([T5l). pQ]) and (p4| immediately give the single 
photon S'-matrix element 



Sp-k = (1 + rk)5kp + TkS-kp , 
with the reflection amplitude 

-iV 

~ 2 J sink {sk ~ n) + iT ' 



(30) 



(31) 



Here, we have utilized DOS p{uj) — 2/-\/4J^ — (uj — wq)^ 
and ReT,{uj) = for the scattered photon in the T-type 
CRA. The result (|30p is the same as that obtained in 
Ref.4. 

In the following, we study the single photon bound 
states in this system by considering the poles of 5'-matrix 
(see Fig. [SK). The poles of ^-matrix are the roots of the 
equation 



0. 



(32) 



FIG. 5: (Color online) (a) The gray shade circle denotes the 
scattering matrix element of the single photon, and it is also 
the Green function of the two level atom. The poles of the 
atom Green function determine the energies of the bound 
states, (b) The band structure of T-type CRA is shown in 
the right panel. 



Because the energies of bound states are not inside the 
photonic energy band of the CRA, i.e., \Eb — uio\ > 2 J, 
the self-energy is 



F sigi-jEb - ^o) 
^{Eb-Lo^f-AJ^ 



Then Eq. ^ becomes 



Eb-n 



T sigri^Eb - ujq) 
^{Eb-ujaf-U^ 



= 0. 



(33) 



(34) 



The above equation possesses two solutions {Elf' and 

-E^"') for Eb'. one solution E^^ is below the bottom of 

energy band, the other solution i?^"' is above the top of 
energy band. The structure of the energy spectrum is 
shown in Fig[5)D. The two bound states are 



\B±)^[Y^rlj±{x,)a\+a+]\Q) , 



(35) 



with the wave-functions^^ in the spatial representation 



(-)i"iy 



V 



E, 



(u) 



id) 



where 



K±{E) 



^ E - ujQ ^^ 



2J 



1 ± 



UJQ- E 
2J 



(36) 



(37) 



It is obviously that In n±{E) < for the both cases 
E — e'^'^ and E = Ef; the wave functions exponen- 
tially decay as |a;| increases as shown in Eq. (1361) . which 
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FIG. 6: (Color online) The wave-functions of two photonic 
bound state in the discrete coordinates. The norms of both 
wave-functions exponentially decay as \x\ increases. 



insure that the photon is indeed localized to form the 
bound states. The wave-functions 'il'±i^i) &re illustrated 
in Fig. ini 

By short summary for this section, we verify the ex- 
istence of photon bound states in the T-type CRA by 
considering the poles of the S'-matrix element. In the 
next section, using the T-type CRA to simulate the T- 
type waveguide, we study the multi-photon transport in 
the T-type waveguide. 



V. SCATTERING MATRIX FOR PHOTONS IN 
THE T-TYPE WAVEGUIDE 

In this section, we focus on the photon transport in the 
T-type waveguide. For the T-type waveguide, S'-matrix 
element is firstly calculated out by the LSZ reduction. 
Secondly, by using the obtained S'-matrix, we give the 
out-state for arbitrary incident state of photons. Finally, 
by analyzing the photon out-state in the spatial represen- 
tation, we can obtain the quantum statistical properties 
of the scattered photons, such as the photon bunching 
and anti-bunching. 

As shown in Sec. II, the waveguide is simulated by 
CRA in the limits k ~ ±7r/2 and ujq — ttJ, in which 
the dispersion relation of the photon is = i;„ |/c| = |A:|. 



Here, we let the group velocity Vg 
Then the Hamiltonian 



1 for convenient. 



Ht = n\e){e\+J2\k\aiak 

k 

V 



-—Y^{ala -HH.c), 



describing the T-type waveguide becomes 

Tji^) _ rre , tto 

with two parts 

fe>0 



(38) 

(39) 
(40) 



k>0 

+-^E(«L^"+H.c.) 

fc>0 



(41) 



where the operators 



Ofc.e 



ak,o 



-^(ofc -I- a-k) 



71 ^"'^ 



(42) 



describe annihilations of the e-photon and o-photon. 
Here, e-photon depicts the photon with even parity in 
the momentum space, i.e., a-k,e — afe.e, and o-photon 
depicts the photon with odd parity in the momentum 
space, i.e., a^k.o — ^o,k,o- In Eq- (|4T|) . the effective cou- 
pling constant is = \/2V . For the o-photon, the is 

diagonalized in the bases {aj; ^ |0)}, so we only need to 
find out the S'-matrix for the the scattered e-photon. In 
the following, we consider the S-matrix for the e-photon 
according to the approach. 



A. Single photon scattering 

For the single e-photon case, Eqs. and give 
the single photon Green's function 



G{p- k) = [Go(fc) 



^-^Glk)]5,k, (43) 



where a — Q, — iTt/'^ and Vt — . Here, we used 
i?eE(a;) = and p(w) = 1. Together with Eq. (fTS]) . 
Eq. gives the single e-photon T-matrix 



p;k 



Spk 



-iVj 



k 



(44) 



Next, we achieve the single e-photon S-matrix element 
Sp-k — tk^pk by considering Eq. (fT3)) . where the trans- 
mission coefficient is 



tk 



(45) 



This result (j45|) accords with that of Refs.— based on 
the Lippmann-Schwinger formalism. 

B. Two photon scattering 

For the two e-photon case, Eqs. (|19p and give the 
two e-photon Green's function as follows: 



G 



pi,P2;fci,fe2 



. 2r| G'o(fci)Go(fc2)Go(pi)Go(p2) 
'(27r)3 (^-a)(c.l-«) 



(wi - a){uj2 - a) 



(46) 



By taking the photon frequency uji and uj[ on shell, we 
obtain the T-matrix 



iT, 



2a) 



TT {p2 - a){ki - a) 
~ (pi " a){k2 - a) ' 



(47) 



From the above equation, the two e-photon S-matrix el- 
ement 



5. 



PiP2;feifc2 



*^PiP2;*:i,fe2 

+ tkitk2{SpikiSp^k,2 +Sp2kiSp^k2) ■ (48) 



follows Eq. (|14p immediately. If two incident photons are 
prepared in the state |fci,fc2), the wave function^ 



{Xc, X \out) 



2 ^ 

PlP2 



SpiP2;kik2 {Xc,x\pi,p2) 



2tt 



[tkitk2 cos(Afca;) 



^p2 ^i{E-2n+irT)\x\/2 



{E-2n + irTy 



(49) 



of two outgoing photons in the spatial representation is 
obtained in terms of the two e-photon center of mass 
coordinate Xc — {xi + X2)/2 and two photon relative 
coordinate x = xi — X2- Here, the total momentum 
(energy) is E ^ ki + k2 and the relative momentum is 
Afc = (fci — k2)/2. When the photon momenta ki and 
k2 both satisfy the resonance condition ki = k2 — fl, 
the envelop wave-function exponentially decays as 
the relative coordinate x increases. This reflects that the 
outgoing two photons attract with each other effectively 
and form a two photon bound state. If £' — 2f2 is kept 
to zero, the wave- function at x = decreases as |Afc| in- 
creases, which implies the photons repulse against each 
other effectively through interacting with the two- level 
atom. The above results about two photon scattering 
accord with the results reported in Refsi^i^.. 



C. Three photon scattering 

For the three photon case, Eqs. ([15]) and give the 
connected T-matrix as 



3(27r)' 



a=l,2,3 



where the functions F*-"-* , are defined as 



(50) 



]\i5(t^i - kp. 



PQ 



{uj[ - a;i)(a;^ - ujy.){uJi - a) 



(ojg — a) [uJi ^ ijJ2~ w'l — a) ' 



(51) 




I 



'Ww/ /WW- — I— 



(a) 



(h) 



FIG. 7: (Color online) Feynman diagrams for three photon S 
matrix: There exist three kinds of disconnected diagrams (a), 
(b) and (c). 



F 



(2) 



^ (,.,L 



PQ 



(w^ - W2)(^ - W3)(a;^ - a) 



(52) 



and 



^.(3) _ 



51 



PQ 



(^2 - W2)(a;i - UJi){ui'i - a) 



U^H^i-PQ,) 



{uj2 - a){uj2 + W3 - u}'2~ a) 



(53) 



Here, P = {Pi,P2,P3) and Q = {Qi,Q2,Qi) are two 
different permutations of (1,2,3), and i = 1,2,3. The 
three photon S'-matrix element 



'Pi,P2,P3\ki,k2,ks, 



* E E E ^Pj\kiTp^,pfj-k^,ks 
i,j=l,2.,3',,5=£i\,l3=£j 

+ E n ^PQ^■''^r'^ 

PQ 1=1,2,3 

+ *Tpi,P2,P3;fel,fe2,fc3 I 



(54) 



is constructed by summing up the contributions from 
all disconnected Feynman diagrams as shown in Fig. [T] 
Here, 7, (5, A, /3 G {1, 2, 3}. The similar discussions can be 
applicable to deal with the TV- photon scattering process. 

Next, we consider the physical meaning of S'-matrix. 
To this end, we first analysis the T-matrices, which 
are \T2\ = |rpip2;fei,fe2 1 and IT3I = |rpi,p2,p3;fei,fe2,fe3 1 • 
Here, \T2\^ describes the two-photon background fluores- 
cence, which is explicitly discussed in Ref.^'^, and jTaj 
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(a) (b) 

FIG. 8: (Color online) Three-photon background fluorescence 
for the total energy of incident photons E = ki + k2 + ks = 3 
and Ft = 1, where the energy level spacing fl is taken as 
units: (a) the three photon are both on resonance with the 
atom, i.e., fci = ^2 = fea = Q,; (b) the energies of the three 
photons are fci — 0.5, fc2 — 0.3, and ks = 2.2, respectively. 

depicts the three-photon background fluorescence, which 
is shown in Fig. |S1 It is shown that when the three 
photons are all on resonance with the atom, the three- 
photon background fluorescence describing jTaj^ is en- 
hanced largely. 

The out-going state of three photons 

\out) = ^ Sp.^p^^p^-kiMM \Pl,P2,P3) , (55) 

Pl<P2<P3 

is determined by the three photon S'-matrix, and its spa- 
tial representation the wave-function reads as 

{xi,X2,X3 \out) 

6(27r)3/2 ■ ^^0^ 

P1P2P3 

The contour maps of the probability distributions 
\{xi,X2, X3\out)\ are numerically shown in Fig. [S] It 
is illustrated in Fig. [9^ that when the three photon are 
all on resonance with the atom, the scattered photons 
prefer the two photon bound state rather than the three 
photon bound state. That is, if two photons form the 
bound state, it is difficult to form three photon bound 
state, namely, the two bounded photons repulse another 
one effectively. When the three photon are not on res- 
onance, the probability distribution \{xi,X2,X3 \out) |^ is 
shown in Fig.[3]D. It is illustrated in Fig.lHb that it is also 
difficult to realize the three photon bound state. If the 
position of one photon is given, such as X3 = 0, other two 
photons do not always attract or repulse each other, but 
attract each other at some points and repulse each other 
at other points, which are determined by the distance 
between the two photon and another photon localized at 
^3 = 0. It is follows the above discussion that, in the 
conceptual setup of the photon transistor, the two- level 
atom controls the coherent transport behaviors of single 
photon, such as the transmission and reflection. In the 




-1 1 

Xf 

(a) 




-1 1 



(b) 

FIG. 9: (Color online) The probability distribution of three 
photon with one photon at origin, where the total energy of 
incident photons E = ki + k2 + ks — 3, the energy level 
spacing Q is taken as units, and Ft = 1: (a) the three photon 
are both on resonance with the atom, i.e., ki = k2 — kg = f2; 
(b) the energies of the three photons are fci = 0.5, k2 = 1.5, 
and ks — 1, respectively. 



multi-photon transport, the atom can induce the effec- 
tive interaction of photons. We can control the effective 
interaction by adjusting the energy level spacing of the 
atom. Therefore, the coherent manipulations of TLS can 
result in a transitions from the repulsive case to attrac- 
tion of effective photon interactions. 

In addition, we point that the recent references ob- 
tained the same results for two photons transport, but 
our works are different from them: (1) Refs.*'^ only give 
the two photon eigenstates by Bethe-ansatz method, but 
we find a general method^i^, i.e., the scattering Bethe 
ansatz technique (see Appendix), to derive the multi- 
photon eigenstates; (2) though we can obtain the multi- 
photon eigenstates by the subtle scattering Bethe-ansatz 
method, we still need a lot of complicated calculations to 
achieve the S'-matrix by using the Lippmann-Schwinger 
scattering theory. However, the LSZ approach can be 
generalized to study the multi-photon scattering in the 
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two level atom 



CRA-1 




H-Type CRA 



FIG. 10: (Color online) The schematic for the H-type CRA 
is shown in this figure. The red circle denotes the two level 
atom. The blue dots denote the coupled resonators. 



waveguide, such as the three photon transport in the 
waveguide; (3) Except these, the LSZ approach is also 
used to deal with the multi-photon scattering in the more 
complex CRA, such as the H-type CRA. In the next sec- 
tion, by using the H-type CRA to simulate the H-type 
waveguide, we investigate the multi-photon transport in 
the H-type waveguide. 



VI. TWO PHOTON SCATTERING PROCESS IN 
THE H-TYPE WAVEGUIDE 

In this section, we study the two photon scattering 
process in the H-type waveguide. The conventional H- 
type waveguide is simulated by the H-type CRA (Fig.[TO|) 
in the high energy limits. The model Hamiltonian 

Hh = r!|e)(e|+ ^ S^Vsal^a- + h.c.) (57) 

'i,s=l,2 

+ XI i^o^A.s'^Us- Js{al,ai+i^s + h.c.)\, 

i,s=l,2 

of the H-type CRA is defined by the hopping constant 
J a and the creation operator a\ ^ of the i-th single mode 

cavity with frequency Wq in the CRA-s, where Vs is 
the hybridization constant of localized atom-photon in 
the 0-th site of the CRA-s. Here, s denotes the CRA-1 
or the CRA-2 as shown in Fig. [TOl In the A:-space, the 
Hamiltonian (|57l) becomes 

Hh = 0|e)(e|+ ^ 

fc,s=l,2 

+ E (Vaal^a- +h.c.), (58) 



k,s=l,2 



where the dispersion relation of photon is e 
2 Jc cos k. 



In the high energy limits k — !• ±7r/2 and ljq = ttJ^, 

the dispersion relation of the photon is e"^ ^ Vs \k\ with 
the group velocity Vs = 2Js. In the case, the Hamilto- 
nian (|58[) describes the photon transport in the H-type 
waveguide. For convenient, we use the operators 



''k,e 
''k,o 



1 

1 

71 



(afc,s + a-k,s) ■ 



(59) 



to rewrite the Hamiltonian ([55|) as Hh — h'^^ + h'^\ 
where 



H 



H 



2^ Efc "fc,, 
fc>0,s=l,2 



k,o ' 



(60) 



E 



fc>0,s=l,2 



for the o-photon and 

H^H^ = ^^|e)(e| 

+ 4f E (V;agV+;i.c.), (61) 

V ^ fe>0.s=l,2 

for the e-photon with the effective coupling Vg — V^Vg- 
Because Hamiltonian H^^ is diagonalized in the bases 
{a^i^^J |0)}, so we only need to find the S'-matrix for the 
e-photon. 



A. The LSZ reduction for e-photon scattering in 
the H-type waveguide 

For calculating the multi-photon S'-matrix, we only 
consider the Green's function and the 5-matrix for the 
e-photon in the H-type waveguide. In this case, the 2n- 
point photonic Green's function reads 

G^;:;;;:^;^;::::X(*i'--4;ii,---^n) (62) 

The S'-matrix element is the overlap 

out{f\i)i^ = m(/|S|i)i„ , (63) 

of incoming wave = '^ki'^e ■■■'^k"'^e \^) ^-"^^ outgoing 
wave state As shown in the Sec. Ill, the basic 

part of the S-matrix is the T-matrix. The relation 



between the 2n-point T-matrix element 



(H) _ rpji,...,j,i;ii,...,i 



pi,...,p„;fci,...,fe„ 



(64) 



(65) 



11 



and the 2ri-point Green's function 



'-^2n — ^pi,...,p„;fci 



:l . (66) 

is determined by the LSZ reduction approach, where 

i 



4'-^ + zO+ 



(67) 



Utr 



is the Green's function of the free photon and 



G 



Pi,...,p„;fei 



is the Fourier transformation of Eq. 
62|) . Here, the subscript os denotes the on shell limits 



LO ^ S 



(0 



Then, the 2n-point Green's function 



pi,...,p„;fci,...,fc„ 



(-l)"j2ninZ[,7|.,r;r 



is obtained by the generating functional 



dujdk- 



where we have used the matrix 

Af[e,r] = 

and the field variable 



4'^ + »0^ 



(68) 



(69) 



As) 



(70) 



(71) 



Here, the atom decay rate is Fg = X]s ^//"s- Finally, 
together with Eqs. §EI and Eq. §^ gives the 

basic element of the S'-matrix, i.e., T- matrix. Using the 
T-matrix elements, we can find the all S'-matrix elements. 
For convenient, we let vi — V2 — 1 below. 



B. Single and two e-photon 5*- matrices 

We consider the single and the two e-photon S- 
matrices in this subsection. The diagrammatic analysis 
shows that the single e-photon S'-matrix element 



consists of the T-matrix element 



(72) 



(73) 



where k and p are the momenta of the incoming pho- 
ton in the waveguide-i and the outgoing photon in the 
waveguide-j respectively. Together with Eq. (|68| and 
Eq. (|69p. we obtain the single e-photon T-matrix ele- 
ment as 



k-n- 



5fep • 



(74) 



Next, we consider the special case: the incident e-photon 
is prepared in the waveguide- 1. Then, the S-matrix ele- 
ments are S^ll = t^^^^Skp and S^.'^ = t^^^^Skp, for 



k-n + q{vi + v,^)' 



and 



.(21) 



k-n 



We can verify the unitarity of the S-matrix as 



(11) 



(21) 



1 



(75) 



(76) 



(77) 



For the two photon case, the diagrammatic analysis 
and the LSZ reduction approach give the two e-photon 
S-matrix element 



PlP2;kik2 Plfel P2fe2 P2fel Plfe2 



-•T^Jl J2;«l,»2 /70\ 

'"piP2;fei,fe2 ' 



with the two photon T-matrix element 



Pl:P2\kl,k2 



(79) 

where fc^ and pr are the momenta of the incoming pho- 
tons in the waveguide-v and the outgoing photons in the 
waveguide- respectively. Together with Eqs. (1681) and 
(|69p . the two e-photon T-matrix element becomes 



■rpjl,j2;il,l2 _ 

PiP2;fei,fc2 



i%%yoiyo2{ki + fc2 - 2ag) 
7r(p2 - a.H){ki - an) 

^ ^ki+k2.,pi+P2 

{pi - aH){k2 - oih) ' 



(80) 



with an = ^ — iTe/^. 

Finally, we consider the special case: the two inci- 
dent e-photons are prepared in the waveguide- 1 and the 
waveguide-2 respectively, i.e., the incident state \in) is 
|0). In this case, Eqs. dTS]) and ^ give 



"fci,e"fe2,e 



cl,l;l,2 
PlP2\kik2 



the S-matrix elements as 

,(ll),(21),c r , c c N 

ki k2 ["kipiOk2P2 ~t~ "kip2"k2Pi ) 

^ tV2Vi^ (fci +k2~ 2aH) 
TT ip2 ~ aH){ki - an) 

^ ^ki+k2,pi+P2 

[pi ~ - an) ' 



S, 



1,2;1,2 
PlP2\kik2 



,(ll),(22)r c ,(21) ,(21) r r 

'•fci ''k2 "kipiOk2P2 + T^ki ^k2 "fclP2"fc2Pl 

iV^^Vi (fci + k2- 2aH) 



IT (P2- Off) (fci 
^ki + k2,pi+P2 

{pi - aH){k2 - an) 



oih) 



(82) 
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FIG. 11: (Color online) The transmission coefficients in the 
waveguide-1 and the waveguide-2: Vi — 1, V2 = 2. The 
energy level spacing Q, is taken as the unit. 



and 



a2,2;l,2 ,(21),(22).r c , r r n 

PlP2-Mk2 ~ fei fe2 VkxV\'^k'2Vl ^ ^kxpi'^k^vi) 

^ iViVj (fci + fc2 - 2ag) 
TT (P2 - aH)(fci - "h) 

'^'£l+fe2,Pl+I'2 



(pi - 0LH){k2 - an) 



(83) 



Here, tj,^^'' is defined by 



t 



(22) _ fc-» + zA(y^2^^^2) 



(84) 



k^n + i^{v^' + vi) 

C. The state of the outgoing photons 

With the help of the S'-matrix elements obtained in 
the above subsection, we find that for the single incident 
e-photon prepared in the waveguide-1, the out-state of 
the scattered photon is 



(ll)^(l)t , ,(21)^(2)t 



\out)^{tral\'J+t- 



"■k.e 



(85) 



As the functions of the incident momentum fc, the t^^^^ 
and 4^^^ are plotted in Fig. [HI We find that if Vi = 
Vi , the transmission coefficient t^,^^'' equals to zero when 
k = Vl (Fig. [T2|) . This result displays that the outgoing 
e-photon is only emitted from the waveguide-2 when the 
incident e-photon is on resonance with the atom. 

For the case of the two incident e-photons prepared in 
the different waveguides, i.e., the state \in) of the incident 
photons is 



\^n) = aW^af t 



|o) 



(86) 



From Eqs. (|8m83p . we obtain the out-state of the scat- 
tered photons as 



\out) = \out)T^^ + \out)^2 + |oui)]^ 



2 ' 



(87) 



FIG. 12: (Color online) The transmission coefficients in the 
waveguide-1 and the waveguide-2: Vi = V2 = 2. The en- 
ergy level spacing SI of the atom is taken as the unit. If the 
incident photon is on resonance with the atom, the transmis- 
sion coefficient in the waveguide-1 equals zero and that in the 
waveguide-2 equals one. this result displays that the outgoing 
photon only emits from waveguide-2. 



with three parts: (a) the state of the two outgoing e- 
photons both in the waveguide-1 is 



\ou 



f] - \^ ci,i;i,2 „(i)t„(i)t |n\ ■ 

Pl<P2 



(b) the state of the two outgoing e-photons in the differ- 
ent waveguides is 

\out),, = io) ; (89) 



(c) the state of the two outgoing e-photons both in the 
waveguide-2 is 

\out),, = ^ Sl^l!l,^a^la^l |0) . (90) 

Pl<P2 

In the spatial representation, the state of the outgoing 
e-photon 

\out) = J dxidx2e''^''^[gii{xc,x)ai\^lai\^l 
+512(2:0, a:)aWt42)t 

+522(2:e,x)a(2)t42)t]|o) , (91) 

is determined by the wavcfunction in the center of mass 
frame Xc = {xi + X2)/2 with the relative coordinate x — 
Xi —X2' (a) the wavcfunction gii{x) of two photons both 
in the waveguide-1 is obtained as 



511 



i(^) = ^[C^4f cos(A.a;) 



4Y2Y3gi{E/2-aH)\x\ 

' AAl - {E - 2aH)^^ 



(92) 
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(b) the wavefunction (712(2^) of two photons in the differ- 
ent waveguides is obtained as 



512(2;) 



4A^ 



{E - 2aHY 



(93) 



(c) the wavefunction 322(2;) of the two photons both in 
the waveguide-2 is obtained as 



522(2;) 



[4f)4f cos(A,x) 



27r 



{E-2aHf 
where we define E = k2 + ki and A^ = (fci — fc2)/2 



(94) 





(b) 



FIG. 13: (Color online) The correlations of the two photons 
in the center of mass frame: Vi — V2 = 2. The energy level 
spacing of the atom is taken as the unit, (a) The total 
energy of the incident photons E is 2, and the difference 
of two incident energies is zero. In this case, the two photons 
are both on resonance on the atom, (b) The total energy 
of the incident photons E is 2, and the difference Afe of two 
incident energies is 1. 



D. Quantum statistics by second order correlation 
functions 

Finally, we analyze the quantum statistical features 
of the scattered photon by the second order correlation 
functions of photons. The second order correlation func- 
tion of the outgoing e-photon is defined by 

G|f {x^,X2) = {out\ a^^la^^la^la^l, \out) . (95) 

It is straightforward to prove that the second order cor- 
relation is just |(7ij(x)|^ by substituting Eq. (|?T|) into Eq. 
(|95|) . Thus the wavefunction gij{x) displays the quantum 
statistic characters of the scattered photons. Here, the 
second order correlation function |5ij(2;)|^ are plotted in 
Fig. [T3] and Fig. [U for different system parameters. In 
these figures, we take the energy level spacing of the atom 
as unit and the total energy of the incident two e-photons 
equal to two, i.e., two times of the energy level spacing 
of the atom. 

Fig. [13] shows that if t^i = V2 =2, the outgoing two 
e-photons attract with each other through the interac- 
tion with the atom. This displays the obvious bunch- 
ing behavior of photons. In Fig. [T3k . the two inci- 
dent photons are both on resonance with the atom, so 
\gii{x)f_= |522(2;)|^ Fig. [Til shows that when Vi ^ V2, 
such as Vi = 1 and V2 = 2, if the outgoing two e-photons 
are emitted from the same waveguide they attract with 
each other and display the photon bunching behavior. If 
the outgoing two e-photons are emitted from the different 
waveguides, they display the photon bunching behavior 
when the two photons are both on resonance with the 
atom. 

By adjusting the difference AF = V2 — V1 between the 
interactions Vi and V2 , we find that (a) If the two pho- 
tons are both resonance with the atom, the out-going 
photons always display the photon bunching behavior, 
however, when AV increases the photon bunching behav- 
ior becomes vague and in the limit Ay 00, \gii{x)\'^, 




-3 3 -3 3 

X X 

(a) (b) 



FIG. 14: (Color online) The correlations of the two photons 
in the center of mass frame: Vi — 1 and V2 = 2. The energy 
level spacing Q of the atom is taken as the unit, (a) The total 
energy of the incident photons E is 2, and the difference At of 
two incident energies is zero. In this case, the two photons are 
both on resonance on the atom, (b) The total energy of the 
incident photons iJ is 2, and the difference of two incident 
energies is 1.8. 

1522(2;)!^ and |5i2(2;)|^ all tend to constants which implies 
the photon bunching behavior disappears, (b) If the two 
photons are not resonance with the atom and the total 
energy equals to 20, the functions |(7ii(a;)|^, 1522(2;)!^ and 
1512(2;)!^ exhibit the oscillation behaviors. As AV — >■ 00, 
\gii{x)\ , 1322 (2;) I and 1312(2;) I all tend to constants and 
the oscillation behavior disappears. 

To summarize this section, we have used the LSZ re- 
duction approach to investigate the single photon trans- 
mission and the quantum statistic properties of the two 
photon transport in the H-type waveguide. The trans- 
port of multi-photon in the more complex CRA can be 
investigated systematically by the same method. 



VII. SUMMARY 

In conclusion, we have demonstrated that the LSZ 
approach is feasible to be generalized for studying the 
complex nano-structure for multi-photon transport, such 
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as the multi-photon transport problems in the complex 
CRA systems. Concretely, the single-, two- and three- 
photon transports in the T-type waveguide are investi- 
gated by the LSZ approach in details. Some of our re- 
sults accord with the known results obtained by other 
methods for the simple case. Besides, the scattering of 
three photons in the T-type waveguide and the quan- 
tum statistical characters of two photons in the H-type 
waveguide are systematically studied by this approach. 
For the single photon transmission in the T-type CRA, 
we find two bound states: the lower bound state is be- 
low the bottom of the photon energy band, the upper 
bound state is above the top of the photon energy band. 
We also gives the three-photon scattering wave function 
by the LSZ approach. Obviously, the LSZ reduction ap- 
proach can be carried out to investigate the iV-photon 
transport in the more complex architectures constructed 
by the T-type CRA. 

Finally, we emphasize that the LSZ reduction ap- 
proach can be extended to investigate not only the multi- 
photon scattering problem in the quantum optics, but 
also the scattering processes in other fields, such as the 
Kondo scattering in the condensed matter physics and 
the Feshbach resonanc o^^i'^^ in the atomic physics. It 
is more interesting to simulate the scattering processes 
in the condensed matter physics or the atomic physics 
by the complex CRA architectures. By making use of 
the LSZ reduction approach, we can study the photon 
scattering in these artificial CRA architectures to under- 
stand the realistic scattering processes as well as electron 
transporl^STi^ in the condensed matter physics. 
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Appendix: Multi-photon eigenstates in the T-type 
CRA 



In the appendix, we utilize the scattering Bethe 
ansata2ii2^ to obtain the multi-photon eigenstates. For 
convenient, we start with the Hamiltonian (1411) . In the 
real space, the Hamiltonian (PT|) is rewritten as H^p = 
Hp + Ha + Hint , where the Hamiltonian of photons is 



Here, 



Hp = —i dxip'' ix)dx'(p{x) 



(A.l) 



and the Hamiltonian of two- level atom is Ha ~ il\e) (e| 
The interaction Hamiltonian Hint is 



H,nt =V dxS{x){ij{x)a^ + h.c.) 



(A.2) 



(A.3) 



is the Fourier transformation of ak^e- The single photon 
eigenstates are constructed by 



\^p) = [ J dxfp{x)i;\x) + epa+] |0) . (A.4) 
The Schrodinger equation H^ = p \ (f>p) gives 

\p) = [ dxe'P^alix) |0) . (A.5) 



where we define an operator 

alix) = [e{-x) + e'^''0{x)]^j^{x) + d{x)epa+ , (A.6) 

which is the single photon creation operator ip^{x) when 
V tends to zero. However, when V is not zero, Q.jj{x) 
neither satisfies the bosonic commutation relation nor the 
fermionic commutation relation. Here, Cp and fp{x) are 



V 



" P-Q + tT/2 ' 
fp{x) = e'P^[e{~x)+e'''-0{x)], (A 

where T = and the phase shift 

p~n-ir/2 



(A.7) 



p-n + ir/2 



(A.9) 



is the same as tk 

For studying the multi-photon eigenstates, we use the 
scattering Bethe-ansatz to assume the A^-photon eigcn- 
state 

N 

\<^>}^ = ^Ap [Dxy^^'-^^^l[alix,)\0} . 

P i=l 

(A.IO) 

Here, [Dx] denotes Hj^i ^(^i+i ~ Xi)dxi and P — 
{Pi, P2, Pn) is a permutation of (1, 2, A^). The 
Schrodinger equation gives the relation Ap/Ap> = 
gj$(Pj ,Pj+i)^ where the phase shift is 



iV 



(A.ll) 



where P = (Pi, P2, •-, -P, , ^^j+i , -, ^w) and P' = 
(Pi, P2, Pj+iiPji Pn)- As an example, we give the 
explicit expression of the three photon eigenstate 



N=3 



/ 9(x3 - X2)9{x2 - xi) dxi 
•' i 

[1 + 5*12 + 5*23 + >5'i2S'l3 

+ 523^13 + 5i25i3^23]e*''^"'+'''""'+'^'"" 

x«l(^i)«l (^2)43(^3)10) , (A.12) 
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where and Pijf{..., ki, kj, ...) = /(..., kj, k^, ...). Here, / is 

any given function. 



kj kj 
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